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ABSTRACT 

The infrared-to-X-ray (IRX) flux ratio traces the relative importance of dust cooling to gas cooling 
in astrophysical plasma such as supernova remnants (SNRs). We derive IRX ratios of SNRs in the 
LMC using Spitzer and Chandra SNR survey data and compare them with those of Galactic SNRs. 
IRX ratios of all the SNRs in the sample are found to be moderately greater than unity, indicating 
that dust grains are a more efficient coolant than gas although gas cooling may not be negligible. 
The IRX ratios of the LMC SNRs are systematically lower than those of the Galactic SNRs. As 
both dust cooling and gas cooling pertain to the properties of the interstellar medium, the lower IRX 
ratios of the LMC SNRs may reflect the characteristics of the LMC, and the lower dust-to-gas ratio 
(a quater of the Galactic value) is likely to be the most significant factor. The observed IRX ratios 
are compared with theoretical predictions that yield IRX ratios an order of magnitude larger. This 
discrepancy may originate from the dearth of dust in the remnants due to either the local variation 
of the dust abundance in the preshock medium with respect to the canonical abundance or the dust 
destruction in the postshock medium. The non-equilibrium ionization cooling of hot gas, in particular 
for young SNRs, may also cause the discrepancy. Finally, we discuss implications for the dominant 
cooling mechanism of SNRs in low-metallicity galaxies. 

Subject headings: galaxies: individual (LMC) — galaxies: ISM — dust, extinction 


1. INTRODUCTION 

A dusty medium is known to be cooled either by an 
infrared (IR) cooling mechanism via gas-grain collisions 
referred to as dust cooling or by a cooling mechanism via 
atomic processes referred to as gas cooling. In hot plasma 
(> 10 6 K), the cooling pro cess is expected to b e domi¬ 
nated by dust cooling (e.g.. lOstriker fc Silk 19731). Since 
the canonical temperature of the swept-up gas in a super¬ 
nova remnant (SNR) exceeds 10 6 K (e.g.. iHughes et all 
119981 ). SNRs provide an ideal laboratory to examine dust 
cooling in the interstellar medium (ISM). 

A measure of the relative strength, the IR-to-X-ray 
(IRX) flux ratio, is often used as an indicator of the 
dominant cooling mec hanism in dusty plasma (iDwekl 
ll987tlDwek et alJll987l) . Using IR data obtained by the 
Infrared Astronomical Satellite (IRAS), IRX ratios of 
nine Galactic SN Rs and four SNRs in t he LMC have 
been m easured bv lDwek et al.1 (1 19871 ) and iGraham et all 
(11987D . respectively, and both studies found that IRX 
ratios of the SNRs are significantly greater than unity, 
supporting dust grains being the dominant coolant. 

As dust and gas cooling pertains to the properties of 
the ISM such as dust-to-gas ratio (DGR) or metallic- 
ity, IRX ratios in different galaxies or in different en¬ 
vironments may reveal a distinct tendency. If an SNR 
or its progenitor has encountered inhomogeneities in the 
ISM or local variations of dust abundance in the ambi¬ 
ent medium, e.g., SNRs in tera cting with nearby molec¬ 
ular clouds dDwek et al.l 119871) , IRX ratios may devi¬ 
ate significantly from the theoretical prediction based on 
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global values. Also, the variation of IRX ratios could 
be a p robe of t he ongoing destruction of dust (e.g., 
iDwek et al.ll2008h . 

The LMC is the nearest gas-rich galaxy to our Gal axy, 
and its metallicity is a half of the Galactic value ('e.g. JPeil 
[1992D . In this context, a comparison of the IRX ratios 
of the LMC SNRs to those of the Galactic SNRs may 
confirm environmental effects on the dominant cooling 
mechanism in the hot plasma. However, no systematic 
comparison of IRX ratios between SNRs in the two galax¬ 
ies has been made so far, probably due to the insufficient 
IR and/or X-ray data on the SNRs. 

In the last decade, the Spitzer Space Telescope and 
the AKARI infrared satellite have enabled large-scale 
IR surveys to be performed towar d the Galaxy as well 
as the Magellani c Clouds (e.g., iMeixner et al.l 120061 : 
IQnaka et al.fl2007h . Using these data, statistical stu dies 
abou t the IR emission of SNRs in the Galaxy (e.g ., Lei 


2005yReach et_al.1 l2006t iPinheiro Goncalves et al 
Jeong et al.ll2012[ ) and the LMG (e.g.. iSeok et al 


2011 


2008 


2013) have been carried out. Moreover, the Chandra X- 


ray Observatory provides superb X-ray data on a large 
number of SNRs such as the Chandra SNR Catalog^. 
The IR and X-ray surveys on SNRs allow us to carry out 
a systematic comparison of IRX ratios and to investigate 
environmental effects on this ratio. 

In this paper, we present IRX ratios of LMC SNRs us¬ 
ing Spitzer IR fluxes and compare them with those of 
Galactic SNRs. Theoretical IRX ratios for the LMC and 
the Galaxy are calculated according to dust and gas cool¬ 
ing functions taking into account relevant DGRs. The 
IRX ratios of the LMC SNRs show a subtle deviation 
from those of the Galactic SNRs, and for some cases, the 
observed IRX ratios are significantly lower than the the¬ 
oretical prediction. We discuss the the possible origin of 
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the observed trend and its implications. 

2. IRX RATIO OF A DUSTY PLASMA 
2.1. Measurement of IRX ratio 

We define the observational IRX ratio as 

IRX ratio = Fi R /F x , (1) 

where Fir, is the total IR flux and F x is the Chandra 
X-ray flux in 0.3-2.1 keV. We adopt t he Spitzer 24 and 
70 /jm fluxes of the LMC SNRs from ISeok et al.l (120131 
hereafter SK013) and X-ray fluxes from the Chan dra 
SNR catalog. Among the 47 SNRs listed in lSKQ13l 29 
SNRs show IR emission at 24 /im, and all but three SNRs 
are also detected at 70 /jm. In addition, X-ray properties 
(gas tem perature and density) of 28 SNRs are listed in 
ISKQ131 and 24 out of the 28 SNRs have IR flux mea¬ 
surements. However, note that we include 19 out of the 
24 SNRs in this study for the following reasons: IR and 
X-ray emission of two SNRs (N157B and N158A) are 
dominated by their pulsar wind nebulae, so dust and gas 
cooling mechanisms are not applicable. For two other 
SNRs (DEM L72 and DE M L256), although th eir X-ray 
properties are presented in lKlimek et al. (2010; their to¬ 
tal X-ray fluxes are not given in the paper. Finally, in 
the case of the SNR in N159, it is inevitable that its IR 
and X-ray emissions are contaminated by a nearby H II 
region and LMC X-l. 

Among the 19 SNRs, we estimate Fir of 16 of them 
by fitting their spectral energy distribution (SED) with 
a single-te mperature modified blackbody (Equation (1) 
in lSKOlOl) . Uncertainties of Fir mostly result from the 
uncertainties of the photometry, which is ~ 15%-20% of 
the total flux. In the case of the three SNRs (0509-67.5, 
0519-69.0, and SN 1987A) for which 70 /jm fluxes are not 
given in 1SKQ131 their Fir are derived on the basis of the 
spectrum taken by the Infrared Spectrograph (IRS) on 
board Spitzer (iDwek et al.l 1 2 01 (j I Williams et al.l 12011L 
see more in Table [T]). 

For all the SNRs but two (DEM L205 and DEM 
L241), we adopt X-ray fluxes from the Chandra SNR 
catalog (see Table [TJ. Adopting the soft-band (0.3-2.1 
keV) fluxes from the Chandra catalog, we expect to ex¬ 
clude a hard X-ray component mainly dominated by syn¬ 
chrotron emission. In fact, the soft-band fluxes are nearly 
the same as the wide-band (0.3-10.0 keV) fluxes except 
for (young) Crab-like SNRs where synchrotron emission 
dominates their SEDs (such as N157B and N158A that 
we likewise exclude above). Although there are three 
Crab-like SNRs (0453-68.5, SNR in N206, and DEM 
L241) in our sample, all are mature enough for their SNR 
shocks to have encountered the ambient medium and to 
produce thermal emission contributing most of the X- 
ray emission. For a few (young) SNRs such as N103B 
or SN 1987A, however, the X-ray spectra indicate multi¬ 
ple components, possibly including a contribution from 
SN ejecta, which need to be distinguished (see section 
14.411 . Uncertainties of the X-ray fluxes mostly arise when 
correcting the interstellar absorption in general, so that 
the X-ray fluxes of the LMC SNRs should be relatively 
secure compared with those of the Galactic SNRs. We 
suppose that the X-r ay fluxes are uncertain by a factor o f 
a few, at most (e.g.. ISeok et al.1120081 : IDwek et al.ll2010t) 
and assume 30% of the X-ray flux as its uncertainty en 


bloc. 

To compare LMC and Galactic SNRs, we accumu¬ 
late IR flu xes o f G alactic SNRs in the literatur e. Re - 
cently, iPinheiro Goncalves et al] (20 IK hereafter IPG fil l 
detected 39 out of 121 SNRs within the Spitzer MIPS- 
GAL survey at 24 and 70 pm and measured their fluxes. 
For 13 of these 39 SNRs, the total IR fluxes and IRX ra¬ 
tios are also measured. We consider that the 70 pm fluxes 
of two SNRs (Kes 79 and G337.2-0.7) are quite uncertain 
due to their low brightness in the MIPS 70 pm. so they 
are excluded from our study. Moreover, we notice that 
most SNRs with the measurement of IRX ratios in iPGIll 
are known to interact with molecular clouds, which do 
not usually show a spatial resemblance between IR and 
X-ray (see below). To supplement Galactic SNRs for 
this study, eight Galactic SNRs including the historical 
SNRs such as those of Kepler or Tycho, are additionally 
included. In total, the re are 19 Galactic SNRs used in 
this work. Since IPG 111 derive IRX rati os using Chandra 
wide-band fluxes (see Table 6 in lPGllh . we newly derive 
their IRX ratios using the soft-band fluxes for consistency 
of this paper although the difference is negligible. 

Furthermore, as the spatial resemblance between IR 
and X-ray morphologies (IRX resemblance, hereafter) 
implies a physical association between dust and hot 
plasma, it is necessary to scrutinize IRX ratios taking 
the IRX resemblance into account. We classify all SNRs 
in our sample into three categories (definite, partial, and 
lack of resemblance) based on visual inspection of the 
Spitzer 24 pm and Chandra soft-band images and use 
the classification for further analysis. 

In summary, we estimate IRX ratios of the 19 LMC 
SNRs and 19 Galactic SNRs as listed in Table [Q The 
total IR and X-ray fluxes, plasma temperatures, and the 
IRX re semblances are given together. Among the LMC 
SNRs, iGraham et al.l ( 1987 ) previously measured IRX 
ratios of N49, N49B, and N63A using the IRAS data 
(~ 12, 4, and 12, respectively). While our new measure¬ 
ment for N49 (N49: 11.1 ± 2.9) is consistent with the 
previously derived ratio, the IRX ratios of the other two 
SNRs (N49B: 1.51 ± 0.39 and N63A: 2.96 ± 0.76) differ 
somewhat from them. This deviation is most likely due 
to the contribution from the FIR emission (see section 

H3). 

2.2. Dust and gas cooling function 

The dust cooling function ( Ah , erg cm 3 s -1 ) in thermal 
equilibrium was computed bv IDwek ( 198 7). and we adopt 
the simpl e analytic expression o f the collisional heating 
rate from lDwek fc Arendt ( 199 2). For a single species of 
dust grains with a given radius (a in pm), the cooling 
function is thus written as 

Ad(T) ~ 2.15 x 10 -21 Tl ,2 h{a , T) (2) 

where p is the mean atomic weight of the gas (p = 1.4), 
Z& is the DGR, T e is the gas temperature in units 
of 10 6 K, p is the material density of the dust grain, 
and h(a, T ) is the heating efficiency. The heating ef¬ 
ficiency is the fraction of electron kinetic energy de¬ 
posited to dust, which is a function of electron energy 
and grai n size ( and species). It was extensively calcu¬ 
lated bv lDwekl (|1987f ). and we adopt the approximate for- 
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TABLE 1 

Observed X-ray and IR fluxes of SNRs 


SNR 

-PlR a 

+x b 

IRX Ratio 

T e c IR-X-ray 

Resemblance d 

(io- 

11 erg cm 2 s - 1 ) (10 

— 11 _0 

erg cm 

s- 1 ) 

(10 6 K) 


LMC Sources 

0453-68.5 

3.84 + 0.69 

1.8 

2.12 + 0.57 

3.4 


DEM L71 

3.85 + 0.66 

2.6 

1.50 + 0.39 

7.5 

D 

N23 

3.57 + 0.56 

2.2 

1.60 + 0.41 

6.5 

D 

N103B 

9.17 + 1.36 

5.6 

1.63 + 0.41 

11.6 

D 

0509-67.5 

0.85 ± 0.43 

0.59 

1.43 ± 0.78 

23.2 

D 

0519-69.0 

2.10 + 1.27 

3.57 

0.58 + 0.38 

17.4 

D 

N132D 

26.3 + 3.7 

17.0 

1.55 + 0.38 

9.5 

D 

N49B 

12.1 ± 2.0 

8.0 

1.51 + 0.39 

4.2 

D 

N49 

79.2 ± 13.5 

7.2 

11.1 ± 2.90 

5.1 

P 

DEM L205 

5.81 + 0.77 

0.048 

121 ± 29.2 

2.9 

L 

SNR in N206 

1.11 + 0.19 

0.44 

2.52 + 0.66 

4.6 

L 

DEM L238 

0.57 + 0.10 

0.072 

7.97 + 2.11 

4.1 

L 

SN 1987A® 

0.31-1.55 

0.12-0.58 

2.4-3.8 (2.5) 

5.0 

D 

N63A 

76.1 ± 12.1 

26 

2.96 + 0.76 

9.0 

P 

DEM L241 

6.57 + 1.46 

0.13 

49.8 ± 14.9 

4.5 

L 

DEM L249 

5.07 + 0.93 

0.046 

110 + 30.0 

8.7 

L 

DEM L316B 

3.41 + 0.66 

0.17 

20.7 + 5.73 

7.5 

L 

DEM L316A 

2.84 + 0.53 

0.12 

24.0 ± 6.55 

16.2 

L 

0548-70.4 

2.04 + 0.34 

0.27 

7.55 + 1.96 

7.5 

P 

Galactic Sources 

Kepler 

131 (1) 

64.4 

2.0 

27.6 (8) 

D 

Gil.2-0.3 

1100 (2) 

395 

2.8 

6.7 (9) 

D 

G15.9+0.2 

530 (2) 

95 

5.6 

10.4 (10) 

D 

Kes 73 

990 (2) 

175 

5.6 

5.8 (11) 

D 

Kes 75 

270 (2) 

19.5 

13 

8.0 (12) 

D 

3C 391 

2000 (2) 

61.7 

32 

6.5 (13) 

L 

W44 

200 (2) 

8.31 

240 

5.8 (8) 

L 

3C 396 

300 (2) 

18 

16 

8.1 (14) 

L 

3C 397 

810 (2) 

129 

6.2 

2.9 (15) 

P 

W49B 

3000 (2) 

587 

5.1 

2.8 (8) 

L 

Cygnus Loop 

7900 (3) 

1400 

5.6 

2.7 (8) 

D 

Cas A 

2570 (4) 

1870 

1.3 

17.4 (16) 

D 

Tycho 

668 (1) 

176 

3.9 

9.9 (8) 

D 

IC 443 

8100 (3) 

1930 

4.2 

2.2 (8) 

L 

Puppis A 

9000 (5) 

3800 

6.9 

7.0 (17) 

D 

G292.0+1.8 

287 (6) 

204 

1.4 

7.7 (18) 

D 

Kes 17 

455 (7) 

5.4 

84 

8.8 (19) 

L 

RCW 103 

2700 (2) 

1700 

1.6 

3.5 (20) 

D 

G349.7+0.2 

1400 (2) 

30.3 

42 

8.8 (21) 

P 

References. — References for Fir and T e 

of the Galactic SNRs. (1) Gomez et al. (20121); (2) PG11 

and references 

therein: (3) Dwek et al. (1 ( 

)87l); (4) Sibthorpe 

et al. (120103); 

(5) Arendt et al. (2010): (6) Lee et al. (2009); (7) Lee et al. 

(20111): (8) Kawasaki et al. 

(20053): (9) Roberts et al. (20033): (10)iRevnolds et al. (2006): (ll)IKumar et al. (2014); (12) 

Helfand et al.l (2003|): 

(13) 

Chen et al. (2004!!); 

(141 Su et al. (20111; (15) Safi-Harb et al. (20051: (16) Hwang Laming 


(2012); (17) Hwang ct al. (2008); (18) Park et al. (2004); (19) Gclfand ct al. (2013); (20) Gotthclf ct al. (1997); (21) 
ILazendic et all (120051 


Note. — 

a Total IR flux (-Fir) of the 19 LMC SNRs and the 19 Galactic SNRs. For all LMC SNRs except three, Fir is derived 
using the Spitzer 24 and 70 /mi fluxes from (SK013) in this work. Fir of 0509-67.5 and 05 1 9-69.0 are derived from 
their dust mass and temperature measured based on their IRS spectra (Williams et al. 2011; Borkowski et al. 2006). 
Uncertainties of their fluxes incorporate the dust temperature range in (Borkowski ct al. (20061). For SN 1987A, see 
below. For the Galactic SNRs, Fir is taken from the references given in parentheses. 

b Soft band (0.3 -2.1 keV) X-ray flux from the Chandra SNR Catalog. Uncertainties are assumed to be 30% en bloc 
(see Section I2.1D , which are taken into account for estimating the uncertainties of IRXs. For two of the LMC SNRs, 
XMM-Newton fluxes are adopt ed; Assu ming the d istance of 50 kpc, we derive the X-ray flux of DEM L205 based on its 
X-ray luminosity measured bv IMaggi et al.l (120121) . For DEM L241, the X-ray flux of the entire S NR is derived by th e 
summation of the X-ray fluxes for its head and tail structures and a point source inside the SNR (|Bamba et al.|[2006l) . 
For four Galactic SNRs, their X-ray fluxes are not included in the Chandra Catalog, or the flux integrated over the 
limited area is on ly listed. In these cases, we too k the X-r ay fluxes from the lit erature; Cygnus Loop and Puppis A: 
IDwek et al.l (119871) , IC 443: IKawasaki et al.l (120051) , Kes 17: IGelfand et al.l (120131) 

c Gas temperature (T e ) in 10 6 K. For the LMC SNRs, T e is taken from (ISKQ13I) . For the Galactic SNRs, T e is adopted 
from the references given in parentheses. 

d Spatial resemblance between IR and X-ray morphologies (D: definite, P: partial, L: lack of resemblance). 
e Due to the youth of SN 19 87A, IR and X-ray fluxes evolve with time. We adopt ranges of the measured fluxes and 
the corresponding IRXs from [Dwek et ahl (|20101 ) between day 6067 to 7983, and the representative value of the IRX is 
2.5. 
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mula given bv iSmith et al.l ( 19961 see the i r equ ation 6). 
We take Z& = 0.0062 from Zubko et al.l (|2004 BARE- 
GR-S model) for th e Galaxy and Z& = 0.0017 from 
iMeixner et all (l2013i l for the LMC. We consider graphite 
as a common dust species for t h e Ga laxy and the LMC 
(p = 2.24 g cm” 3 , iZubko et al.ll2004 1 although taking a 
different dust species suc h as silicate does not alter our 
results much (see Section [4711) . 

Figure Q] shows the dust cooling functions for a single¬ 
sized dust populations with radii of 0.05, 0.1, and 0.2 
(im (equation [2j . A<j for the Galaxy in Figure Q] shows 
good agr eemen t with A d from more elaborate calcula¬ 
tions by iDwek et~aLl (l2Q08h . In addition, gas cooling 
functions are shown in the figure for comparison. For the 
gas cooling, we compute the cooling curve in collisional 
ionization equilibrium (CIE) by using the CHIANTI code 
1V7.0: lLandi et, aTl 1 201 2ft . The elemental abundances of 
the G alaxy and the LMC are taken from lAsolund et,~aI1 
(|2009ft and iRussell fe Dopital (I1992T) . respectively. Note 
that we calculate not only the total gas cooling but also 
the soft X-ray (0.3-2.1 keV) cooling curve, Ax, which 
we adopt for more direct comparison with the Chandra 
X-ray fluxes. 

Using the cooling function of a dusty plasma via gas- 
grain collisions, A d (T), and of the gas via atomic pro¬ 
cesses, Ax(T), the theoretical ratio of the IR cooling to 
X-ray cooling, hereafter referred to as the IRX t h ratio, is 
defined as 

IRX th ratio = A d (T)/A x (T). (3) 

It is worthwhile to note that there are several caveats 
in comparing the ratio of theoretical cooling rates to the 
observed IRX ratios. One would result from the fact 
that heavy elements are repeatedly included when we 
calculate the gas and dust cooling rates. Behind the SNR 
shock wave, dust grains are destroyed and some fraction 
of heavy elements might be returned to the gas phase. 
So in SNRs, the DGR would be lower while the metal 
abundance in the gas phase would be higher than those 
in the general ISM. We assume the non-depleted LMC 
(solar) abundance for the gas cooling, considering that 
the X-ray emitting gas will have an abundance close to 
it. If some fraction of heavy elements are locked into 
dust, the gas cooling curve will be lower than those in 
Figure [lj We will further discuss other caveats in Section 

ECU 

3. RESULTS 

Figure [5] shows the IRX ratios of the 19 LMC SNRs 
and 19 Galactic SNRs as a function of plasma temper¬ 
ature and the distributions of the IRX ratios. The IRX 
resemblance is also denoted for each SNR in the figure: 
definite, partial, and lack of IRX resemblance are shown 
respectively by solid, open with cross, and open symbols. 
Based on A d and Ax in Figure [TJ the IRX t h ratio for the 
LMC and the Galaxy is shown in Figure [2] 

First of all, the SNRs in both galaxies (except 0519- 
69.0) show that IRX ratios greater than unity, which 
supports the idea that dust cooling is the primary cooling 
mechanism in the dusty plasma. However, it is found 
that more than half of the LMC SNRs have IRX ratios 
just between 1 and ~ 3 (see also Table [TJ . Moreover, IRX 
ratios considerably greater than unity are mostly derived 
from those without a definite IRX resemblance. This is 



Fig. 1.— Cooling function of a dusty plasma via gas-grain col¬ 
lisions for a single-sized dust population with radii of 0.05 (dotted 
line), 0.10 (dot-dashed line), or 0.20 fim (dashed line) from Equa¬ 
tion ©. For comparison, cooling curves via atomic processes in 
the soft X-ray (0.3-2.1 keV) regime (Ax: thick solid line), and to¬ 
tal gas cooling curves (A g : thin solid line) are overlaid. Red lines 
represent cooling functions for the LMC while blue lines are for the 
Galactic conditions. 

somewhat surprising in the sense that IDwek et al.1 (|1987ft 
report that most of the Galactic SNRs have IRX ratios 
much larger than unity. In the case of 19 Galactic SNRs 
shown here, one third of them have IRX ratios below 
~ 3, seven have ratios between 3 and 10, and the rest of 
them have ratios even larger than 10. 

It is noticeable that the LMC SNRs with a lack of the 
IRX resemblance appear to have generally higher IRX 
ratios (median ratio of those with the partial or lack of 
IRX resemblance: 21) than those with definite IRX re¬ 
semblance (median: 1.5). A similar trend can be found 
from the Galactic SNRs, too (median ratio: 3.8 for those 
with definite IRX resemblance and 24 for the rest). These 
SNRs with high IRX ratios (or without definite IRX re¬ 
semblance) are mostly known to interact with molecu¬ 
lar c louds (e.g., 3C 391, W44, or IC 443, l.liang et, ahl 
120101 ) and/or are so-called “mixed-morphol ogy” SNR s 
with centrally brightened X-ray emission dRho fc Petrd 
I1998H . For SNRs interacting with molecular clouds, the 
IR emission can be enhanced because radiation may play 
an important role in heating dust grains compared with 
collision; this eventually results in high IRX ratios. This 
will be further explored for the Galactic SNRs in our 
future paper. 

Among the SNRs with definite IRX resemblance, the 
IRX ratios of the LMC SNRs seem to be systematically 
lower than those of the Galactic SNRs (recall the median 
IRX ratios of 1.5 and 3.8, respectively). Such a difference 
is, in fact, consistent with the trend seen in the theoreti¬ 
cal expectation that the IRX t h ratio for the LMC is lower 
by a factor of ~ 2-3 than that for the Galaxy although we 
need to be cautious about the direct comparison between 
the measured IRX and IRX t h ratios. The comparison 
between the observed IRX and IRX t h ratios is relevant 
only to those with definite IRX resemblance because the 
rest of them might have contributions from radiatively 
heated dust. We compare only those IRX ratios with 
the theoretical prediction. As shown in Figure [2j all the 
LMC SNRs with definite IRX resemblance have much 
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Fig. 2.— Left: IRX ratios of LMC (circles) and Galactic (squares) SNRs as a function of plasma temperature (T e ). Based on their IR 
and X-ray morphologies, the SNRs are classified into three groups with definite, partial, and lack of resemblance designated by solid, open 
with cross, and open symbols, respectively. Owing to lack of space, two Galactic SNRs (Puppis A and G292.0+1.8) are not named in the 
figure. Theoretical IRX ratios for single-sized dust populations with radi i of 0.05 (dotted lines) and 0.20 pm (dashed lines) are overlaid 
for the LMC (thick red) and Galactic (thin blue) conditions (see section [2.21 for the explanation). Right: a histogram showing the IRX 
distributions. The IRX distributions of the LMC and the Galactic SNRs are denoted (thick red and thin blue dashed lines, respectively) 
together with the total IRX distribution (solid line). 


lower IRX ratios than IRX t h ratios by up to an order of 
magnitude. Such a large deviation has been previously 
r epor ted for a few Galactic SNRs, too (e.g.. iDwek et al.l 
ll987t ). and this could be att ribut ed to the low DGR or 
dust destruction (see Section ITTil) . 

4. DISCUSSION 

4.1. IRX Ratios of SNRs in the LMC and in the Galaxy 

While theoretical calculations predict that dust cooling 
plays a primary role in the cooling of a dusty plasma, the 
relative importance of dust cooling (A<j) to gas cooling 
(Ax) is expected to depend on environmental properties 
such as DGR, metallicity, and grain compositions. Our 
sample shows that the measured IRX ratios of the LMC 
SNRs are systematically lower than those of the Galactic 
SNRs, and the theoretical IRX ratios can explain this 
observed trend to some extent. The IRX t h ratio of the 
LMC is lower than that of the Galaxy (Figure O, which 
is due to the low DGR of the LMC. As shown in Figure 
Q] while Ax at T > 10 6 K is reduced by < 50% due 
to the half-metallicity of the LMC (mainly Fe, Mg, and 
O), A<j of the LMC is four times lower than that of the 
Galaxy due to the 1/4-DGR of the LMC. The measured 
DGR of the LMC might have some variation depe nding 
on m ethods or observational data (e.g., iMeixner et al.1 
l2010t ), but nevertheless the lower IRX ratios of the LMC 
SNRs relative to those of the Galactic SNRs seem to 
be most significantly attributed to the low DGR of the 
LMC. 

We also examine other possibilities that may system¬ 


atically influence IRX ratios. Grain properties such as 
composition or size distribution may be different between 
the LMC and the G alax y though we assumed the same 
properties for both. iPeil ( 1992 ) examined the extinction 
curves of the LMC and the Galaxy and suggested that 
dust in the LMC is composed predominantly of silicate 
while the abundances of silicate a nd carbonaceous g rains 
are comparable in the Galaxy. IMeixner et al.l (20101 ) also 
suggested that the standard properties for Galactic dust 
are not appropriate and that amorphous carbon is re¬ 
quired for the LMC instead of graphite. However, the 
densities of different grain species do not significantly al¬ 
ter the dust cooling function (Equation ©), which can 
be a factor of two at most (e.g., p = 1.81-1.87 and 
3.5 g cm~ 3 for amorphous carbon and silicate, respec¬ 
tively; iZubko et al.ll2004l) . In addition, the details of the 
grain size distribution in the LMC are not well known, 
but the relative abundance of very small grains to large 
grains in the LMC is rather similar to the Galactic value 
((Bernard et al.l[2008D . This may indicate that there is no 
significant difference between the LMC and the Galaxy 
in terms of the grain size distribution. Therefore, the dif¬ 
ferent grain composition and size distribution could be 
the subsidiary factors to the low IRX ratios of the LMC 
SNRs. 

4.2. IRX ratio and IRX resemblance 

The SNRs without definite IRX resemblance are found 
to have high IRX ratios. This may imply that there is 
some (or more) contribution to the IR emission origi- 













































6 


Seok et al. 


nating from radiatively heated dust or atomic/molecular 
emission lines rather than collisionally heated dust. For 
the radiatively heated dust, there are two aspects. In 
the postshock gas, radiative heati ng might domina t e dust 
heating. Using the Spitzer data, lAndersen et al. ( I201lh 
modeled the dust emission of 14 Galactic SNRs and com¬ 
pared the relative importance of collisional heating to ra¬ 
diative heating as a function of grain size (see Figure 5 
in their paper). In most cases, it is found that radiative 
heating is dominant for grains with a size of 4.0-110 nm. 
This is not surprising given the fact that these SNRs 
show signs of interaction with molecular clouds, which 
induces radiative shocks to produce t he strong UV ra¬ 
diation from the postshock gas (e.g., iDraine fe Me Keel 
119931 : IKqoI I2014H . Another aspect is that dust grains 
in th e presh oc k gas can be heated by shock radiation 
(e.g., iLee et all 12011 1. If the preshock medium is dense 
gas such as a large molecular cloud, then a significant 
fraction of the shock radiation will be absorbed and con¬ 
verted into IR radiation whereas the X-ray emission from 
shocked gas is only a fraction of shock radiation. This 
can naturally explain the high IRX, and in this case we 
do not expect to see any spatial resemblance between IR 
and X-ray emission. 

Additionally, ionic and/or molecular line emission 
could contribute to Fir. For example, in N49, one of 
the SNRs that show strong emission lines in the IR spec¬ 
trum, the emission lines can contribute ~ 38% of the 
MIPS 24 pm flux, and the [O I] 63 qm line pr ovides 11% 
of the MIPS 70 pm flux (lOtsuka et all 1201(1 and refer¬ 
ences therein). However, since this fraction is not enough 
to attribute the high IRX solely to the contribution from 
the line emission, the line emission could be a relative ly 
minor cooling source as concluded by lAndersen et al.1 

(MU). 

4.3. Comparison to the theoretical IRX 

Figure [2] clearly shows that the LMC SNRs with a defi¬ 
nite IRX resemblance, have a significant deviation of the 
IRX ratios from the IRX t h ratios. iDwek et~aI1 (11987T ) 
found a similar discrepancy from a few Galactic SNRs 
such as Kepler or SN 1006. Both these SNRs are located 
far above the Galactic plane and have much lower IRX 
ratios than the predicted values. This may suggest that 
the severe dearth of dust in the ambient medium aro und 
the SN Rs results in the discrep ancy. [Borkows ki et al.1 
(|2006H and I Williams et all (|2006fl calculate dust destruc¬ 
tion in four Type la and four core-collapse LMC SNRs, 
respectively, which overlap with our sample. They in¬ 
fer the total dust masses swept up by the blast wave 
before sputtering, and these masses are generally several 
times lower than the typical value from the average DGR 
of the LMC. Interestingly, 0453-68.5 (and 0548-70.4), 
which has a high DGR of 9.8 x 10~ 4 (and 7.5 x 10 -4 ), 
comparable to the canonical value of the LMC, shows a 
relatively high IRX ratio whereas the rest of them that 
have lower DGRs show much lower IRX ratios (and con¬ 
sequently a more significant deviation between IRX and 
IRX t h ratios). This implies that the low IRX ratios may 
arise from the local fluctuation of DGR in the preshock 
medium of the SNR. 

Another major reason for the low IRX ratios might be 
because dust grains are destroyed behind the shock front 
and therefore the DGR decreases with time. Dust grains 


swept up by SNR shocks are destroyed by thermal and 
non-thermal sputtering with ions. At gas temperatures 
above ~ 10 6 K and shock speeds > 500 km s _1 , the sput¬ 
tering rate is approximately constant, and the character¬ 
istic lifetimes of carbonaceous (C) and silicate (Mg 2 Si04) 
grains with radius a in pm are a/\da/dt\ « 1 x 10 6 and 
2 x 10 5 a/nn years, respectively, where »h is th e d ensity 
of H nuclei of hot plasm a (see Figure 2 in Noza wa et al.1 
120061 : IDwek et al.1 120081. Therefore, for old SNRs, we 
expect theoretical IRX ratios lower than the curves in 
Figure [2j For young SNRs, there is yet another fac¬ 
tor. Since it takes time for the swept-up ions to reach 
the CIE, the X-ray cooling rate of young S NRs is much 
higher than that in th e CIE model (e.g., IDwek et al.l 
119871 : ISmith et al1ll996H . Therefore, for both young and 
old SNRs, we expect them below the equilibrium IRX t h 
ratio in Figure [2] It is anticipated that theoretical IRX 
curves taking into account the dust destruction and the 
non-equilibrium ionization cooling effects will provide a 
better understanding of the observed IRX ratios. 

4.4. Uncertainties of IRX ratios 

While we carry out the systematic comparison of IRX 
ratios in the two galaxies, there are some uncertainties 
in measuring IRX ratios that we can improve if better 
knowledge of X-ray emission and more comprehensive IR 
analysis of SNRs are achieved. The X-ray emission in an 
SNR may consist of different physical components such 
as thermal emission from swept-up ISM or synchrotron 
emission. In most cases, the Chandra X-ray spectra 
can be explained by the swept-up ISM with a single¬ 
temperature component, which is generally soft (< 1 
keV). We have searched available literature to check the 
X-ray contribution from SN ejecta or multiple compo¬ 
nents that are unlikely to be associated with the regions 
where dust cooling occurs. For some SNRs (e.g., N103B, 
N132D, SN 1987A, and N49), however, their X-ray spec¬ 
tra definitely require multiple components or a contribu¬ 
tion from ejecta (e.g.jLew is et, al.ll 2003l:jPark et, al.ll2003l : 
iBorkowski et, al.lf2007t IDwek et alJl2008lh In such a case , 

the X-ray flux can be corrected, as IDwek et al.1 (120081 1 
applied a correction factor to disentangle hard and soft 
components in the X-ray emission of SN 1987A, but in 
general it is difficult to separate it out. Thus, due to the 
contribution from SN ejecta, IRX ratios will be lower 
than IRXth ratios, in particular for young SNRs. In¬ 
deed, as shown in Figure O young SNRs with ejecta- 
dominated X-ray fluxes are all located i n the lower part 
of the diagram (e.g., N103B: iLewis et al.1120031 SNR 
0509-67.5 and 0519-69 0: IWilliams et all I2011L Gas A: 
iHwang fc L;)iniimll2012l Kenler: iRevnolds et, al.ll2006lf . 

While we estimate the total IR fluxes (Fir), fluxes at 
wavelengths longer than 70 ^m are not taken into ac¬ 
count. If there is a significant far-IR (FIR) emission from 
cold dust in the SNR, we might underestimate the to¬ 
tal IR flux. For example, iGraham et al.l (I1987H include 
IRAS 100 pm fluxes for their estimation. For the two 
SNRs (N49B and N63A) showing some differences be¬ 
tween the new and previous measurements of IRX ra¬ 
tios (Section 12.11) . their IRAS fluxes at 100 pm are the 
strongest among the fluxes in the four IRAS bands. 
However, since both SNRs probably have considerable 
FIR confusion due to nearby bright sources (e.g., the pho- 
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toionized western lobe embedded in N63A and a bright 
clump near the southern limb of N49B, see Figures 4-5 
in ISKQ131 L their previous IRX ratios might be overesti¬ 
mated. 

Recently, the Herschel Space Observatory has provided 
FIR data with un precedented high sensitivity and spa¬ 
tial resolution. iOtsuka et al.1 (1201011 detect the bright 
FIR emission from N49 using the Herschel images and 
show that a significant emission originates from the cold 
dust component. In spite of its superb capabilities, how¬ 
ever, only two SNRs, N49 and SN 1987A, have been 
clearly detected by H erschel so far (jOtsuka et al.ll20l7il : 
iMatsuura et all 1201 If) . Note that FIR emission of SN 
1987A originates not from the swept-up ISM/CS M but 
from the newly formed dust (llndebetouw et all 120141) 
whereas most LMC SNR s in our sample are rat her ma¬ 
ture SNRs (> 10 3 years). lLakicevic et al.1 (12015H present 
the FIR atlas of the LMC SNRs using the Herschel data 
and report that SNRs are not clearly distinguishable in 
the FIR bands. Moreover, no evidence for large amounts 
of newly formed SN-origin dust in SNRs is found, and 
they argue that cold interstellar dust in SNRs is signifi¬ 
cantly destroyed due to sputtering by SN shocks. In this 
context, it is less likely that the FIR emission has a sig¬ 
nificant influence on IRX ratios in this work although we 
cannot rule out the subtle contribution. 

4.5. Implications for low-metallicity galaxies 

The observational trend that IRX ratios of the LMC 
SNRs are systematically lower than those of the Galactic 
SNRs implies that dust grains in extreme low-nretallicity 
galaxies may not be a dominant coolant of hot plasma 
unlike in normal-metallicity galaxies. If the DGR is 
signi ficant ly reduce d as the metallicity decreases (e.g., 
iLisenfeld fc Ferraralll998t lHirashitalll999tl , the dust cool¬ 
ing would be reduced almost proportionally to the DGR. 
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Meanwhile, the gas cooling is modestly affected in hot 
plasma (> 5 x 10 6 K) in proportion to metallicity. 
For example, the SMC has a metallicity even lowe r 
than the LMC (1/5-1/8 Z@, IRussell fe Dopital 1199211 . 
There are more than 20 SNRs in the SMC identi¬ 
fied by X -ray a nd/or radio observations (e.g., 23 SNRs 
listed in I Badenes et al.l l2010f ) , and the recent XMM- 
Newton survey of the SMC reveals that 20 SNRs show 
X-ray emission (lHaberl et al.l [2012lh Surprisingly, how¬ 
ever, only two SMC SNRs (IE 0102.2-7219 and B0104- 
72.3) have been found to have the IR cou nterparts 
('e.g.- IStanimirovic et al.ll2005l : IKoo et al.ll2007lh The de¬ 
tectability of IR emission in the SMC SNRs is extremely 
low whereas about 60% of the LMC SNRs show asso¬ 
ciated IR emission (jSK013lh Although no systematic 
study of the SMC SNRs has been carried out, this might 
suggest that the dust cooling in the SMC SNRs is less 
efficient, most probably d ue to the low DGR o f the SMC 
(Zd = 1.1 — 8.3 x 10~ 4 . IMeixner et al.l l2013h . Further 
IR observations toward the SMC SNRs using higher spa¬ 
tial resolution and higher sensitivity such as NIRCarn or 
MIRI onboard James Webb Space Telescope will help in 
understanding the cooling mechanisms of the SMC SNRs 
or of hot plasma in the low-metallicity environment in 
general. 
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